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This paper deals with an improved backstepping control strategy for 
sensorless control of a permanent magnet synchronous motor (PMSM) based 
on field-oriented control (FOC), using a backstepping controller to improve 
its performances. However, this control requires the precise knowledge of 
some machine’s variables which could not be available. In electric drives 
control, sensors are generally used as the main devices for feedback 
information. Some practical constraints could affect the system 
performances, due to the lack of measurement material or maintenance 
difficulties caused by the dysfunction or faults of the used sensors such as: 
encoder or resolver sensors of speed-position. In this paper, a sensorless 
control is proposed based on a dynamic backstepping method and an 
extended Kalman filter (EKF) which uses the state space formulation with a 
set of mathematical equations to recursively estimate future observations and 
minimizes the mean square error of the estimated variables (rotor speed 
position and torque) to design controllers for nonlinear systems. The 
proposed control scheme achieves the asymptotically uniformed stability. 
The effectiveness of this method is illustrated by the stabilization and 
tracking numerical examples and the obtained simulation results show the 
effectiveness and the feasibility of the proposed controller using Lyapunov 
approach. 
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1. INTRODUCTION 


Permanent magnet synchronous motor (PMSM) is widely used in industrial applications compared 
to other electric motors such as: induction motor (IM), direct current (DC) motor and switched reluctance 
motor (SRM) based on the required characteristics [1], [2]. Mainly, due to its compact design, high 
efficiency, high torque to inertia ratio, high dynamic response, excellent reliability, great robustness, power 
density, technological maturity, a gap for advancement in controller design, small volume, and reduced 
maintenance [3]-[5]. In other side, the non-linearity of the dynamic model of the PMSM produces a large 
specific control difficulty [6], [7]. Parameters, load torque variations also the coupling between the motor 
speed and the electrical quantities, such as the d-q axis currents, make this system obviously hard to control 
[3], [8], [9]. This motor can be controlled by the conventional proportional integral (PI) classical controller 
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[10], [11]. A large techniques and approaches of PMSM nonlinear control have been extremely studied, 
developed and proposed such as: input-output linearization control [12], robust control [13], sliding mode 
control [14]-[16], back-stepping control [17], [18], fuzzy logic control [19], direct torque control (DTC) [20], 
[21], and linear parameter varying (LPV) [22]. But the problem or difficulties facing the used approaches are: 
the proportional integral derivative (PID) controller cannot guarantee satisfactory performance such as: 
stability and control against disturbances. Howover, the problem of chattring of sliding mode controller 
(SMC) requires control that can switch at an infinite frequency. Thus, during the sliding regime, the 
discontinuities applied to the control can train to chattering phenomena, called reluctance or 
"chattering". These switchings deteriorate the control precision, can prove deterioration of mechanical parts of 
the systems and increase the the temperature in electrical systems (significant loss of energy). However, we 
must keep in mind that in fuzzy logic, it is impossible to predict the performance of a system due to its 
dependence on the well selection of the setting parameters. On the other hand, for the DTC control strategy, 
the oscillations of the various variables are sensitive, compared to the current-voltage measurement sensors, 
and limed by the controlled variables: the torque and stator flux. Recently, the back-stepping presents 
promisingly alternative methods for controlling nonlinear systems by combining the choice of the Lyapunov 
function with the control laws which allows the overall stability of the system [17], [23]. 

Since mechanic position and speed sensors (encoder, resolver or hall sensor) are usually too expensive, 
they increase the cost and complexity, decrease the stability of the system, impose hard maintenance, and are 
sensitive to noise. Thus, sensor-less control is becoming very attactive research and a big challenging task [18]. 
The state observer permits an accurate estimation of the speed and the position in the presence of system noise 
and parameters variations. The sensorless PMSM control has been exauxted in the last decades. Trabelsi et al. 
[18] proposed the control of sensorless PMSM using sliding mode observer, Luenberger [23], back- 
electromotive force (back-EMF) [24], model reference adaptive system method (MRAS) has proved [25], or the 
artificial intelligence. Among those, the extended Kalman filter (EKF) is the most well-known and optimal one, 
especially in the presence of the process and measurement noise [26]—[30] which is well adapted to nonlinear 
systems and such applications. In this paper, an optimal PMSM field-oriented control strategy is proposed, 
using a backstepping controller and an EKF based on a recursive method to estimate the rotor speed position of 
the considered machine and to minimize the mean squared error. 

The contribution in the present work focuses on the combination of the EKF with backstepping 
control in order to estimate the parameters or the measurement of the state machine due to the fort non 
linearity of system parametres and the faced difficulties in the sensors, based on backstepping approach for 
better tracking of the speed by introducing the estimation of both rotor speed-position and electromagnetic 
torque respective. The main advantage of introducing such controller is its skills to cope with some 
drawbacks of the PI controller. In practice, each variable requires two controller gains, which are 
continuously changed, due to some of imprecise phenomena in the system. In contrast of that, and as 
developed later, less control variables are used, which guaranties the stability of the system via the choice of 
only three positive gains. Furthermore, an EKF is used to estimate both the rotor speed and position, as can 
be seen in Figure 1. A sensor-less control is introduced to increase the robustness, to ensures better dynamic, 
high performances, and decreases the cost and complexity. 
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Figure 1. The proposed nonlinear sensorless controller of PMSM 
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To achieve the proposed approach, the rest of the paper is organized as; section 2, is reserved to the 
mathematical model presentation of the considered machine. The design of the backstepping controller is 
presented in section 3. Section 4, is dedicated to the design of the nonlinear sensor-less control using EKF. 
The simulation results demonstrate the effectiveness of the proposed control which are illustrated and 
commented in section 5, followed finaly by a conclusion. As perspective, we plan to optimize the developed 
algorithms by the advanced approaches such as: particle swarm optimization (PSO), wolf gray, genetic 
algorithm (GA), neural network, and implement them experimentaly with embedded control devices such as: 
(STM32F407) with its fourth-generation M4 processor and Infineon type intelligent power module (IPM) for 
power part. 


2. PMSM MAMATHEMATICAL MODEL 
The dynamic nonlinear model of PMSM in the (d-q) reference frame is given as [12], [15]: 


dig Rs. Lq 1 
=-—j iQ V, 
dt La a tg P'a Tia d 
Rs. L r p® 1 
—4 = — Sig + pigd-—t2+—V, (1) 
lq lq lq lq 
dQ 3p. 1 


3P sa f 
aE = 2J iq +3 (la — La)iaig — ;2- FT. 


By choosing: x = [ig ig)? = [x, X2 x3 |" as a state vector and the load torque Ty, as an external 


disturbance, u = [ V aVą JË as control vector, Then we can write the state model of the machine in the form: 


Rs 
— 7 PX3 1 
a j Rs por L i Va sl 
X2 =|Px -7 TT Wee ee gs ł v+ La ay She (2) 
X3 apde sf | LX L Bite! X3 
0 a 0 0 


The dynamic model of the nonlinear system is presented with the following representation (3): 


dx _ 
fe ae ee , with: f(x) = Ax and h(x) = Cx (3) 


y =h) 


3. BACKSTEPPING CONTROL DESIGN 

The recursive back-stepping technique guarantees the asymptotic stability of the system in close 
loop, by choosing an adequate Lyapunov function for stabilizing error variables [16], [18], [23], [31]. 
Furthermore, the controller is very efficient for nonlinear systems. In fact, the PMSM dynamic model is 
highly nonlinear, and such control technique can be applied directly with no linearization. The vector control 
principle is generally introduced to assimilate the considered machine control to the DC motor. The ig current 
is always forced to be zero [24] by the orientation of the all linkage flux in the d axis to obtain the maximum 
torque. The principal objective of the proposed back-stepping controller is to ensure the reference speed 
tracking of the PMSM drive in the presence of external disturbances changing (the torque load) [24]. The 
proposed robust controller is designed in two steps is being as: the first step concerns the speed controller 
design methodology and the second, incorporates an imbricate current controller, via an augmented 
Lyapunov function. 


3.1. Speed control design methodology 
The speed loop tracking error variable and its derivative is given respectively by (4) and (5) 


Cx, = X3re¢ — X3 (4) 
ex, = X3 ef — X3 (5) 
From (2) and (5), the speed error derivative becomes: 


f 
eee 7X3 + i (6) 


ex, = X3ref gi 2J 
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By choosing a positive candidate Lyapunov function: 
1 
v = te, (7) 


Its derivative is given by: 


By aca ab Sts Se ; 3p®f f 1 

Vi = €x, x; = €x; (haret = “yy #2 + j%3 + ET.) (8) 
The back-stepping design method considers the d-q axes currents x, andx, as virtual control elements and 
specify its desired behavior, which are called stabilizing function in the backstepping design terminology is 
being as: 


Xref = 


2 With Ką, is a positive constant 9 
ipa; (fxs +T, +J. Kya: €x) xa 1Sa p (9) 


xref = 


In order to guarantee the reference tracking and according to Lyapunov stability theory, the condition V, < 0 
must be verified. Substituting (9) in (8) the derivative of V4, one gets: 


V, = —K,,€2, Kx, > 0 (10) 
Therefore, the speed error approaches zero and global asymptotic stability is achieved. 
3.2. Backstepping current controller 

To ensure, the asymptotic stability of the system (1), a second step is added, allowing the d-q 


currents to match the set values. The following current errors are defined as: 


Pe = Xref AA 


era = Azp — 1 with x;,,, = 0 (11) 


Their dynamics can be written: 


; : ; R 1 
Oxy = Xie — 1 = Ta — px2X3 — zVa (12) 
ROSS bye 2 Rs pof 1 

Êx, = Ž2 ef — 42 = sos +T, +J. Kyz- €x) + T X2 + px1x3 + 5X3 =z (13) 


The stability of the studied system is analyzed based on Lyapunov function (14): 
1 
V, = 5 (Cx, + eg ee) (14) 


Based on (6), (12), and (13) the derivative of (14) is given by: 


V2 = ey, €x, + €x, Ox, + Ox, 6x, (15) 
V, = —k,,e2, —k,,e2 —ky,e2, + ex, |k Hd p Bs + 
2.55 X3 Cx, x1 Ex, X2 ex, ex, x1 &x, TE L X3X2 
2(kx3J-f) 3p®f 3p®f Vq Rs Df 
Gy | Kx, €x se ey ky, €x,) + F ex — F + 1 X2 + X3X1 + X3 TA (16) 


The expression (15) must be defined semi negative, and consequently the control laws are deduced as: 


Va = ky, Ley, + Rsx; — Lx3x2 (17) 
2(kx3J-f) (3p@ 3p®fL 
Va = a ae ob — Kyg@xg) + rex, + RsXz + Lxgxy +X: + ky, Ley, (18) 


With this choice, the derivatives of (14) become: 
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Vz, = —kgeg — kaea — kgeg < 0 (19) 


This means that the tracking errors will converge asymptotically to zero. 


4. SPEED AND POSITION OBSERVER DESIGN USING EXTENDED KALMAN FILTER 
In this section, the EKF observer is applied to estimate the PMSM rotor speed and position in order 
to achieve the feedback process control performed by the previous back-stepping control strategy. 


4.1. Extended Kalman filter principle 

The EKF is a mathematical tool working in iteration and numerical way, able to reconstruct the 
system states from other measurable physical variables [26]. This estimator is an optimal predictor-corrector 
which it minimizes the estimated error covariance when some presumed conditions are occurred [26]. For the 
speed and the position estimation of PMSM, where parameters variation and measurement noise are present, 
EKF is the ideal one [32], [33]. The Figure 2 illustrates the EKF principle. The EKF reposes on some 
hypothesis especially the noise. Effectively, it supposes that the noises which affect the model are Gaussian 
and white and are decorrelated from estimated states [32]. 


i = f (Xkr-1  Up—1) + Wer 


20 
Yr = hey) + Vk (20) 


Where, xp = [ia, tq, Qrp Ox Tr, 1"; The estimated vector; Yk = [hı h2 ]" The system output; ug: the 
control law; wę: noise process; Vg: noise measurement are respectively stand for errors of the parameters and 
measurement; f (Xk-1 , Ux—1), h(x,): nonlinear functions, where: f (Xk-1 , Uk-1) = Xk-1 + X~-1Ts; Ts is the 


i 
sampling time and h(x) = fal The noise covariance matrixes are defined is being as: 
qk 
Q = cov(w) = cov(w,w") = E(w,w’*) (21) 


R = cov(v) = cov(v,v") = E(v,v") (22) 


PMSM model 


Xe = f (Xp-1 »Ue-1) + Wed 
Yr = h(xy) + VE 


Oh(x;,) 
Ox, 


Of (k-11 »Ux-1) 


Oxf, r 
E Xk-1=Xk|k-1 


A(Rrix—1) 


Figure 2. Simplified scheme of the EKF principle 


The choice of the covariance matrices elements is the critical point in the design of EKF which is 
related to the performance and the convergences of the system. The two matrices Q and R depend on the 
training parameters, the sampling time of the amplitudes of the measurements, and some other secondary 
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factors. Also, the robustness of the EKF is linked to the chosen values of the sampling time, which is 
generally linked to the switching frequency and also linked to the capacity of the control board used. In 
practice, the sampling time is linked to the execution speed like the DSP board [34], [35]. The EKF 
estimation procedure is realized in two steps: prediction and correction (update). 


4.1.1. Prediction time update 

In this phase, a prediction of the state vector is given %,\,-1 at sampling time (k) from the input 
Ux and state vector X,_1x-1 at previous sampling time(k — 1). The first time we go through two points of 
initial estimation: Point one: predefine the state measurement noises covariance matrix Q and R. Point two: 
Initialize the filter error covariance matrix P. After the initialization of the estimated matrix then, we have 
two steps must be done in this phase: 
Step 1. Project the state ahead: the predicted state estimate is determined by; 


Rea = f (ĉk-1]k-1 »Ue-1) (23) 
Where f (êk-1k-1 ja) = k-1jk-1 + Ts Ëk-1lk-1 (24) 


Step 2. Project the error covariance ahead: the predicted covariance estimate is calculated as; 


Peja = FkPg-1ijk-1Fk + Qk-1 (25) 
_ Of Êk-1lk-1Uk-1) S _ ôh(xp) 
Where Fẹ = = yg a 2 ; Hk = a E 
Xk-1=ĝk|k-1 Xk=ĝfk]k-1 


4.1.2. Correction-observation phase (update) 

The updated states estimate X,\, which are obtained from the predicted estimated states X;),_, by 
adding a correction term K, Yç to the predicted value. So, three steps are defined: 
Step 1. Compute the Kalman gain Kx: defined by the expression; 


Ky = Prje—1 He (HkPkjk-1Hg + Re)" (26) 


Step 2. Update state estimation with measurement: the expression of the updated state estimate at the output 
of the observer is given by; 


Kite = Êkjkr-1 + Ke We (27) 
Where Yk = Vix — Îkjk-1 (28) 


With, y, = h(x,) : actual output vector; and, Jxix—1 = h(êkik-1): predicted output vector. So, the final 
expression of the update satates estimation is: 


Rete = Rena + Keeler — Ahlâk- )] (29) 
Step 3. Update error covariance is defined as: 
Prk =(I- Ky Hg )Pkik-1 (30) 


4.2. PMSM rotor speed and position estimation using extended Kalman filter 
The PMSM model with state and measurement noises is: 


far = fx) +g.ut+w 


(31) 
y=h)+v 
For the estimation using EKF, (31) is discretized directly using Euler approximation proposed in [25]. 
Xk = f (Xk-1 Ups) + Wed 
32 
Re = A(x) + Vk GA 
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R 1 
(1 -Ts 25) Xip- H TsPX3p 1X21 + Ts Var 


R 1 1 
Moreover: f (xp-1,Ug1) = (1 -Ts “s) Xap — TsPX3y_,%1y4 — Ts 7 PP SX +Ts Vara (33) 
f 3p 1 
(1-154) xara + Ts Bopp Ts Tui 
Nays + TspX3,_, 
xy 1000 0 
And h(x) = By Hy = lo 100 3l GY) 
R 
1—-Ts > TspX3, TsPX2, 0 
R R 1 
1- is 1- Ts = —Tspx1, — Ts 7p Py 0 
So: Fy = (35) 
0 Ty = A i= T L 1 
0 0 Tsp 0 
0.002 0 0 0 0.02 0 0 0 
0 0.002 0 0 0 0.02 0 0 
Qk = ; Rk = (36) 
0 0 0.002 0 0 0 0.02 0 
0 0 0 0.002 0 0 0 0.02 


5. SIMULATION RESULTS, DISCUSSION AND INTERPRETATION 

The simulation results are carried out and perfomed in Matlab-Simulink to validate and verify the 
effectiveness of the proposed sensor-less nonlinear controller with regards to the following tasks: possible speed 
change and load torque disturbance. The used bloc scheme is described in Figure 1. The nominal PMSM 
parameters are shown in Table | and the values of the gains are: ko = 700; kg = 10000; ką = 10000. 


Table 1. The PMS motor parameters 
Stator resistor R,=1,40. 
d and q stator inductances Lg = Lg = L = 0.0058 H 
PMSM pole pairs number p=3 


Permanant flux Pp = 0,1546 Wb 
Inertia J = 0,00176 kg. m? 
Friction f = 0.000388 N.m 
Vdc 400 volts 


The used speed profile refence is shown in Figure 3 which contains four intervales: 

- During the first-time interval [from 0 s to 1.5 s], The set speed changes linearly with time, until it matches 
a low constant value of 50 rad/second. In this case we notice that: the back-stepping controller arrives to 
track the evolution of the reference speed even when the load torque is introduced and the introduced 
EKF permits the estimated rotor speed to track the actual motor speed, presented in Figures 3 and 4. 
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Figure 3. Dynamic speed trajectory tracking 
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Figure 4. The speed tracking error between reel rotation speed and rotation speed estimated 


During the second time window [2 s to 3 s], the proposed controller-observer based strategy in high speed 
is teted where the set motor speed is kept at 400 rad/s. Figure 3, shows clearly the reliability of the used 
strategy as well as its strength and rigidity in all conditions, both for load torque rejection speed regulation. 
The controller arrives to reject the effect the load torque changing (see Figure 3) in a very short period of 
no more than 0.006 seconds, representing the response time. Also, in the steady state, the difference 
between the set and the actual speed does not exceed 2 rad/sec, which represents 0.49%. This is in fact 
enough to say that this strategy is consistent and robust. 
To test of the proposed strategy in severe condition and exacltly during the third time interval [from 3.5 s 
to 4 s], a zero-set speed is chosen. As can be remarked, both the back-stepping controller continues to track 
properly the reference velocity, with no notable tracking error. Furthermore, the estimator continues to play 
the role properly, where the estimated speed matches the measured velocity adequately. This means, that 
the EKF can be used in all the range of the chosen reference speed profile. 
Finally, during the last interval [4.5 s to 5 s], a reverse speed operation is planned. As can be seen, notable 
performances are obtained, where the estimated velocity matches very well the set speed changing in both 
dynamic and steady state. 

Figure 5, illustrates the actual and the estimated rotor position results for the PMSM drive. The error 


between the actual and estimated position of the rotor is due to the use of voltages and currents instead of 
reference voltages and currents. Through the graph Figure 6, representing the three-phase current we find that 
there are two observations, first a direct proportional for current amplitude. In the sense that whenever the 
torque changes, there is a change rapidly as load torque varies. The second concerns the rotation velocity and 
the current frequency is inversely proportional. For the observation mentioned above, it is accurately 
explained by the two curves of the current axis q and torque as shown in the two Figure 7 and Figure 8. 
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Figure 5. Dynamic position trajectory tracking 
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Figure 6. Three-phases Ia, Ib and Ic (A) stator current axes tracking 


20 


TEO S 
= —Iq_est 
3 
= 
= 
g 
© 
OS 
= 
= 
£ s5 
a 
2-10 
£ 
AS | | | | | | | l 
o 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
Time (seconds) 
Figure 7. The iq-axis stator current 
— Torque machine 
——Torque_estimate with EKF 
Torque load _ 
Ẹ 
z 
S 
S 
o 
E 


o 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
Time (seconds) 


Figure 8. The torque tracking 


In the present scenario, the load torque is applied from 0.7 s to 1.3 s and between [2.2 s - 2.9 s], with 
a magnitude of 12 Nm and 15 Nm. The Figure 3 shows very well the efficiency of the estimator. Although in 
the presence of external perturbation, the EKF tracks the actual value of the speed with a very small 
estimated error as shown on Figure 4 with the measured value of rotation speed. The estimated rotor speed 
follows the measured one during all changes in speed profile. According to Figure 9, the effectiveness of the 
proposed nonlinear controller is proved for tracking a reference speed. The igcurrent is forced to zero, which 
matches with the vector control technique idea, and Figure 9 illustrates that. 

In Figure 10, we notice that this error graph evolution over the time between the measured current 
(Id_mes) and the estimated (Id_est) is very small does not exceed the range [—2 * 1077,2 * 1079] it is almost 
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equal to the zero value. Now this proves that, the estimation algorithm the FKE is very reliable and robust 
concerning the parameters of the system such as the load torque or others. Figure 10 shows the measured 
electromagnetic torque evolution over the time span. As can be seen, the motor torque tracks adequately the 
evolution of the load torque, with adequate sign, concerning the speed direction operation. The locus of the 
stator flux Pa-Pg is shown in Figure 11. It can be seen, that in the stator reference frame, the curve is 
elliptic, with no chattering phenomenon remarks, which proves that the both stators flux ®,-®g are 
sinusoidal. 


Id (A) 
—Id_est (A) 


The current Id & Id est(Amper) 


2 2.5 3 
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Figure 9. The Id and Id_est -axis stator current tracking 
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Figure 10. The current tracking error between Id and Id_est (A) 


o 
Flux 2, 


Figure 11. The locus of the stator flux Pa-Pg 


5.1. Test of the robustness of sensorless control with echelon profile of speed 

Figures 12 and 13 show the PMSM speed responses while expediting up the rotor speed from 0 to 
300 rad/s and decreasing up of rotor speed at stopping then changing the rotor direction to -200. The 
reference speed is changed in a period of 1 s with a sequence [50 rad/s], [100 rad/s], [200 rad/s], [300 rad/s], 
and [-200 rad/s]. Figure 12 indicates that the optimal EKF estimated rotor speed is closely tracked against the 
reference speed and overlaps the rotor speed sensor. As the rotor speed increases, the external load becomes 
higher. Therefore, more torque is needed in the engine. The current id is almost zero and has a few small 
pulses when shifting gears. Figure 13 is more clarified with the zoom option at each step. 

Table 2 and Figure 13 illustrate the performance of the system when changing the disturbance of the 
external load and the reference speed. When the machine is started increasing the rotor speed in the sequence 
of [50 rad/s], [100 rad/s], [200 rad/s], [300 rad/s], and [-200 rad/s]. The rotor speed reaches 300 rad/s after 3 s 
(nominal speed). At the instant t=1.25 s we apply an external load of value 5 N.m, the speed of the rotor will 
drop to 99.8 rad/s and stabilize again at the setpoint of 100 rad/s at the instant t=1.252 s. Then the speed 
reduction or the overshoot equals 0.2 rad/s and the recovery time equals 0.002 s. 
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Figure 12. Dynamic speed trajectory tracking robustness 
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Figure 13. Zoom speed trajectory tracking robustness 


The engine speed is reduced to 299.6 rad/s at time t=2.25 s and stabilizes again at the value 
200 rad/s at time t=2.252 s when the load torque is increased to 10 Nm. The decreasing value or the 
overshoot equals 0.004 rad/s and the recovery time is 0.002 s. Even the electromechanical torque also is well 
following the load torque as shown in Figure 14 with small overshoots. We conclude from the obtained 
results that the operation of the engine is always running in a stable manner with this approach sensorless 
backstepping control algorithm using EKF estimator under dynamic load conditions. Table 3 represents the 
nomenclature of the used parameters in mathematical formulas. 


Table 2. Performances of the controller and of sensorless EKF 


Time St Perfo trolle MEER Performances Parameters of 
windows ep COMARES CORTO ER paee sensorless EKF controller 
Maximum Steady Steady Steady 
Plages Settling overshoot Mp state state error state error 
(Sec) Response load torquë time (s) While applying speed (Q_mes- (T_mes- ko ka ka 
disturbing torque (rad/s) Q est) T_est ) 
[0-1] No load 0.0027 0% 49.995 0.02% 0.025% 700 10000 10000 
[1-2] Full load to 5 N.m 0.003 0.2% 99.98 0.03% 0.025% 700 10000 10000 
[2-3] Increased load to 10 0.004 0.15% 199.96 0.03% 0.025% 700 10000 10000 
N.m 
[3-4] Ful load with 10 N.m 0.009 0% 299.9 0.0375% 0.025% 700 10000 10000 
[4-5] Decreased load to 0 0.005 0.07% 0.02 0.01% 0.025% 700 10000 10000 
N.m 
[5-6] Full load with 5 N.m 0.018 0.1% -200.02 0.02% 0.025% 700 10000 10000 
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Table 3. The nomoclature 
B Control matrix Ko, ka andkg Speed, id and iq current controller PI gains 
A System state matrix Cx 41 Oxy) Oxy Id, iq current, and speed error 
C Output matrix Xk The vector state 
w External disturbance vector P covariance matrix State 
6 Rotor position Q covariance matrix noise system 
Q Rotor angular velocity R Covariance matrix noise measurement 
d, q Two-axes synchronous frame quantities Ts Sampling period 
İq» ig-Va, Vą d-q axis current and voltage components k Sampling index 
T, Load torque np Number of pole pairs 0, K Initial value and Kalman filter gain matrix 
igrefs laref__d-q axis reference currents k|k-1 Predicted estimate 


6. CONCLUSION 

In this paper, the nonlinear back-stepping control strategy with EKF based on a recursive method to 
estimate the rotor speed position of the considered machine and to minimize the mean squared error is 
developed. The proposed control approach shows the dynamic speed behavior estimation taking account the 
parameters variation and load perturbation. The performance of the proposed controller has been investigated 
in numerical simulations using MATLAB environment and the presented technique shows a better dynamic 
speed response over the set changing speed, where the introduced controller arrives to reject the effect of the 
load torque in a finite time and the developed speed sensorless was successfully designed with the used EKF 
observer. The obtained results show that, the stochastic estimator arrives to synthesize the rotor speed in all 
conditions with small estimation error. In the future studies, it is planned to design an electronic interface 
where more parameters can be added and changed. Also, similar numerous developments can be done by 
using embedded system. It is also planned to try the training phase of the developed algorithm on more 
power systems and experimental part with the introduction of the advanced control and optimization 
techniques. 
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